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Abstract 

The Wigner function is argued to be the only natural phase space function evolving 
classically under quadratic Hamiltonians with time-dependent bilinear part. This is used to 
understand graphically how certain quadratic time-dependent Hamiltonians induce squeezing 
of quantum states. The Wigner representation is also used to generalize Ehrenfest's theorem 
to the quantum uncertainties. This makes it possible to deduce features of the quantum 
evolution, such as squeezing, from the classical evolution, whatever the Hamiltonian. 


1 Introduction 

The Wigner function [1] can be used to get a visual understanding of why certain time- 
dependent Hamiltonians induce squeezing of quantum states. We first address the question: Why the 
Wigner function, rather than, say, the Husimi function [2]? It will be showed that although other phase 
space functions may evolve classically under certain specific quadratic Hamiltonians, the Wigner 
function is the only one to do so under quadratic Hamiltonians having a time-dependent bilinear part. 
We consider, as an example, squeezing by a periodically modulated harmonic oscillator. We then 
discuss a generalization of Ehrenfest's theorem applying to the quantum uncertainties. This allows to 
deduce aspects of the quantum evolution, such as squeezing, from the classical evolution even in the 
case of arbitrary Hamiltonians. 

Phase space variables (position and momentum) are denoted by q and p. A caret is used to 
identify operators. Thus, [q, p] = i, where we take h = 1. The time variable is denoted by t. 

2 Quadratic Hamiltonians and phase space functions 

Ehrenfest's theorem expresses the time derivatives of the expectations (q) and (p) in a way 
formally similar to Hamilton's equations. In the case of a Hamiltonian H = lp 2 /m + V(q), it reads 

£-<q> = (p)/m, ±<p) = ( F (q)> ( 1 ) 
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where F(q) = -9V/dq is the force. If the potential V(q) is quadratic, hence F(q) linear (or, if the 
particle is in a state sufficiently localized on the scale of non-harmonic variation of V), then (F(q)) = 
F((q)), and (1) become identical to Hamilton's equations, implying that (q) and (p) follow classical 
trajectories in phase space. This stays true in the case of a general quadratic Hamiltonian 

H(t) = H^t) + H 2 (t) (2) 

where Hj(t) is linear, and H 2 (t) is bilinear in q and p(a,p, a, b, c scalar functions of time): 

Hi(t) = Hi(q,p,t) = a(t)q + p(t)p (3) 

H 2 (t) = H 2 (q,p,t) = a(t)q2 + b(t)fop+pq) + c(t)p 2 (4) 

Can one push this further, and associate with state vectors ly), or state operators p, a phase space 
function f(q,p,t) whose quantum evolution is classical for such Hamiltonians? That is, such that 

|j-f(q,p,t) = {H, f} PB = (^^- ^-^fCq.p.t) or f(q,p,t) = f(qt,p t ,0) (5) 

where (qt,pt) is the point which classically evolves into (q,p) in the time interval (0,t), under the Hamil- 
tonian H(q,p,t). It is natural to ask that f(q,p) be linear in either ly) or p, so that it has the form 

f(q,p; y) = <<()qply> or f(q,p; p) = Tr{Aqpp} (6) 

where l<t>qp> and Aqp are kets and operators parametrized by phase space. 

Introduce the unitary time-evolution operator U(t,t'), defined by 

Jj-U(t,t') = -iH(t)U(t,t'), |rU(t,t') = iU(t,t')H(t'), U(t,t’)U(t',t) = U(t,t) = 1 (7) 

and similarly Ui and U 2 corresponding to Hi and H 2 . Let us first assure classical evolution under 
linear time-independent Hamiltonians Hi = aq+Pp: Here, Ui(t,0) = e'^Q+Pp), and one gets classical 
evolution, f(q,p,t) = f(q-pt, p+at, 0), up to a phase, iff 

I0qp) = Dqpl4*)> Aqp = DqpADqp (8) 

where 

Dqp = e'Pq-'QP (9) 

are phase space displacement operators, and 1<J>) and A are some fiducial ket and operator. We now 
note the relations ([ ] + denotes anticommutators, { ) PB Poisson brackets): 
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(10a) 


[Hi, Dqp] Hi(q,p)Dqp, Dqp]+ - [Hi, Dqp} ra = 

i[H2, Dqp] = {H 2 , Dqpjpg (10b) 

where H u = Hi 2 (q,p). Using these results, and refering to (5), one finds that 

U(0,t) Dqp = Dq iPt U 2 (0,t) e>X(q.P.O (1 i) 

where U 2 is the time evolutor for H 2 , and %(q,p,t) is a phase determined by the equation 

T t x - {H(t),x)p B = |Hi(t), x(q.P.O) = 0 (12) 

One verifies (11) by verifying that both sides of U(0,t)DqpU 2 (t,0) = Dq^e* satisfy the same 
differential equation. We then get, from (6) and (1 1): 

fy(q»p; V(0) = <4)IDqptU(t,0)lv> = e-i%(q.P.O f^Cq^p,; y(0)), = U 2 (0,t)l<» (13) 

fA(q,p; p(0) = TrfDqpADqp 1 U(t,0)pU(0,t)) = f^fqt.Pt; p(0)), A t = U 2 (0,t)AU 2 (t,0) (14) 

One sees that the function evolves classically iff l(J>t> is stationary, i.e., if H 2 is time-independent, 
and l<J>) — IE n ) is an eigenket of it. For instance, if H 2 is a harmonic oscillator, and l<{>) = IEo) is the 
ground state, then !<})qp) are coherent states, and 


f <t>(q>p; V( 0 ) = e’ l * (q,p,t)+itEo f<j)(q t ,pt; V( 0 )) 


(15) 


is the Husimi function [2], which is well known to evolve classically under that H 2 of which l<j>) is 
an eigenket. Clearly, no function fy linear in state vectors can possibly evolve classically if H 2 is 
time-dependent. Consider now (14): Again, f A evolves classically iff A t is stationary: For instance, if 
H 2 is time-independent, and A = IE n XE n *l or A = g(H2). then the evolution is classical - but only for 
that specific time-independent H 2 . Is there an operator A (apart from the unit operator) such that 


U2(0,t) A U 2 (t,0) = e* 1 ^ 1 ) A 4# [H 2 (t), A] = <p(t) A 


( <p(t)=3tp/at) (i6) 


(tp a real phase) for time-dependent H 2 (t)? Yes, the parity operator fi, since qf[ = -flq, pit = -ftp 
imply [H 2 , fl]=0. Setting A = it in (14) yields the Wigner function [1,3] 


f w(q,p, p) — K ^Tr{Dqpfl Dqp p], 


f w(q,p; p(t)) = f w (qt,Pt; p(0)) 


(17) 


The Wigner function is the only phase space function which evolves classically under (any) time- 
dependent H 2 (t). One can see this as follows: Represent any A by its Weyl symbol A w (q,p), where 
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the Weyl symbol of an operator A is defined as [3,4] 

A\v(Q>p) = 2Tr{Dqpfl Dqp A} (18) 

We want At, hence its Weyl symbol, to be stationary for any H 2 (t). Now, just as the Weyl symbol 
(17) of p evolves classically under H 2 (t), so does that of At. The only possible way for (A t ) w (q,p) 
to be stationary under classical evolution with different H 2 's (e.g., with orbits which are ellipses or 
hyperbolas of different eccentricities) is that it be be concentrated at the origin: We must thus have 
A w (q,p) = (const)5(q)8(p); this implies [5b] that A = (consOfl. 

3 Squeezing by a Periodically Modulated Harmonic Oscillator 

The classical evolution of the Wigner function under quadratic H(t) is very useful for under- 
standing the quantum in terms of the classical evolution. As an example, consider a harmonic oscillator 

H = Im-Jp 2 + Imco 2 ^ 2 = lto(ap 2 + q 2 /a), a = (moo)- 1 (19) 

If a=l (in suitable units), the classical orbits H(q,p) = constant are circles in the phase plane (q,p). If 
a * 1, the orbits are ellipses (Fig.l), the ratio of the q semi-axis to the p semi-axis being a. We will 
now let a alternate between two values, y and y 1 (where y > 1), at every quarter of a period 2 ti/co, 
while keeping co fixed (i.e., only m changes). Let a = y < 1 for the first quarter period: During that 
time interval, a point initially on the positive q axis (beginning of trajectory 1 on Fig.2) moves to the 
negative p axis, while receding away from the origin by a factor y. Then let a = y > 1 during the next 
quarter period: The point moves to the negative q axis, receding away from the origin by another factor 
y. And so on. One here has parametric amplification. On the other hand, points initially on the p axis 
(trajectory 2) close in on the origin. Thus, classically, the phase plane gets squeezed into the rotating q 
axis, by a factor y at each quarter period. Whence a corresponding squeezing of quantum states [5a,6]. 
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4 Ehrenfest Theorem for the Quantum Uncertainties 

The above concerned quadratic Hamiltonians. It will now be indicated that the Wigner-Weyl 
(WW) representation allows to reinterpret Ehrenfest's theorem, and to extend it to quadratic obser- 
vables, hence to the quantum uncertainties. This makes it possible to deduce features of the quantum 
evolution from the classical evolution even in the case of arbitrary Hamiltonians. Let us first recall that 
in the WW representation, expectation values of operators A have the "classical" form [ 1 ,7] 

<A> = Tr { A p ) = Jdqdp A w (q,p) f w (q,p) (20) 

In particular, the quantum expectations of q and p are 

<q) = jdqdp q f w (q,p), <p> = jdqdp p f w (q,p) (21) 

and the uncertainty matrix is 

Cqq = ((q - <q» 2 > = jdqdp (q - <q» 2 f w (q,p) 

Cqp = Cpq = (J(qp+pq) - <q)(p» = jdqdp (q-(q»(p-<p» f w (q,p) (22) 

C PP = <(P * <P» 2 > = jdqdp (p - (p» 2 f w (q,p) 

A quantum state (wave packet) may be roughly represented in phase space by an uncertainty ellipse 
(u - <u» C’ 1 (u - <u» = 1 where u = (g), C = ^ ) (23) 

We also need the result [7] 

[A, B] w (q,p) = i{A w , B w ) PB (q,p) if A or B is quadratic (24) 

that is: The Weyl symbol of the commutator of two operators, one of which is quadratic, is equal to the 
Poisson bracket of the individual Weyl symbols. Let now A be an observable quadratic in q and p. 
We then have, by (20) and (24), for any H(t): 

&<A> = -iTr{A[H, p]} = -iTr([A,H]p} = jdqdp {A w , H w } PB f w (q,p,t) 

= jdqdp A w (q,p){H w , f w } PB (A quadratic, any H) (25) 

where in the last line we performed an integration by parts. Eq.(25) says that the rate of change of the 
expectation of a quadratic observable is classical, in the sense that it is the same as if each point in 
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Wigner phase space instantaneously followed a classical trajectory. This does not mean that (A) 
evolves classically over finite time intervals, because f w (q,p,t) in (25) is the exact quantally evolved 
Wigner function at time t, not one evolved classically during some finite time (unless H is quadratic); 
more specifically, the first time derivative of (A) is classical, but not the higher order derivatives. By 
letting A stands for q or p in (25), one gets Ehrenfest’s theorem: 

|f(q) = Idqdp q(H w , fwlra = Idqdp {q, H w )pBfw(Q>P>0 = (3H/5p) (26a) 

|j-<p> = Idqdp p { H w , f w ) pb = Idqdp { p, H w } ra fw(q,p,t) = (-dH/3q> (26b) 

The last expressions in (26a,b) are the usual statement of Ehrenfest's theorem [equivalent to (1) if H = 
5 P 1 2 3 4 5 6 7 /m + V(q)], giving a formal quantum-classical analogy. The first expressions in (26a, b) tell us 

much more: That the rates of change of (q) and ( p) are classical, but relative to a phase space 

distribution function. Eq.(25) generalizes Ehrenfest's theorem to quadratic observables, and thus to the 
uncertainties: Indeed, according to (21), (22) and (25), the uncertainty ellipse (23) evolves exactly as if 
each point in Wigner phase space instantaneously followed a classical trajectory. For instance, if the 
ellipse gets squeezed classically, during some small (infinitesimal) time interval, then so does it 
quantally. In general, one may expect that if the classical motion during a finite time interval squeezes 
the uncertainty ellipse, then so does the quantum evolution. The latter statement is of course rigorously 
true if H(t) is quadratic, in view of (17); it is also approximately true, in the case of arbitrary H(t), if 
f w (q,p,t) is sufficiently localized on the scale of non-harmonic variation of H(q,p,t), for the evolution 
of f w (q,p,t) is then approximately classical [5b]. 

Let us mention, finally, that Ehrenfest's theorem for quadratic observables can also be written in 
a form corresponding to the last expressions in (26a,b), namely [5c] 

^<A)= <{A W , H w } PB (q,p)w> (A quadratic, any H) (27) 

where the subscript w on the function [A w , H w ] PB (q,p) of the non-commuting operators q and p 
signifies that they are ordered according to Weyl's ordering rule [4]. 
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